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We consider a SU(3) ring-exchange model on a triangular lattice. Unlike the SU(2) case, under perturbation 
expansion of the SU(3) Hubbard model, the three-site ring exchange is present as well as the usual four-site 
ring exchange. Interestingly, the three-site ring exchange differs from the usual two-site and four-site exchanges 
by a minus sign and is ferromagnetic. We first present numerical site-factorized state studies on this model 
which shows three-sublattice order and a ferromagnetic phase. We further study the model using slave-fermion 
mean field in which we rewrite the exchange operators in terms of three flavors of fermions. We find the main 
competing trial states are the trimer state (triangular plaquette state) and the gapless U(l) spin liquid states with 
parton Fermi surfaces which include both the uniform zero-flux spin liquid state, and the uniform 7r-flux spin 
liquid state. Furthermore, we find there are possible pairing instabilities of the zero-flux (Fermi surface) spin 
liquid state toward a /-wave gapless (nodal) spin liquid state and the 7r-flux (Fermi surface) spin liquid state 
toward an interesting exotic s-wave "gapless" spin liquid state with two flavors of fermions paired up while one 
flavor of fermions remains gapless. 



I. INTRODUCTION 

In recent years, the cold atomic systems have become a 
powerful tool to realize strongly-correlated systems, includ- 
ing N-flavor Hubbard model on different lattices. Such sys- 
tems consisting of several flavors of interacting fermions can 
be realized as different hyperfine states of alkali atoms or 
nuclear spin states of ytterbium 10-01 or alkaline-earth atoms 
Jll-01. A model Hamiltonian to describe such systems is the 
N-flavor fermionic Hubbard model JUS Hi 
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where a, j3 run over the different flavors, (jk) runs over pairs 
of nearest neighbors on the lattice, and j runs over all lattice 
sites. 

The most common case is when N = 2, which corresponds 
to the usual spin- 1/2 Hubbard model. If we focus on the half- 
filling condition, i.e., when each site is occupied by exactly 
one fermion, the system undergoes metal-to-Mott insulator 
phase transition for sufficiently large repulsion U. In exper- 
iments on cold atoms, the transition to a Mott insulator has 
recently been observed. [l^] In the large U regime, it is 
generally accepted that the ground state is well-captured by 
the usual anti-ferromagnetic(AF) Heisenberg model and the 
charge degrees of freedom are completely frozen out. How- 
ever, in the regime where U is not sufficiently large com- 
pared with the hopping strength t, the ground state is not 
well-understood. It is possible that the strong charge fluctu- 
ations play important role for stabilizing "gapless" spin liquid 
phases in this regime. lfTTl[l2Tl 

For the more general case with N > 2, lfl3l[l4ll if we focus 
on certain fillings, the Mott insulating states will also emerge. 
In this case, the spin order is not understood even in the large 
U limit in which we can only focus on the Heisenberg-like 
(two-site exchange) Hamiltonian 



(jk) 



(2) 



where Pjk is so-called two-site exchange operator, which per- 
mutes the fermions between two nearest-neighbor sites as 
Pjk\ot, p) = a), where the a, j3 represent the spin states 
at sites j and k. For N = 3, there has been numerical evidence 
on such SU(3) Heisenberg model on a triangular lattice sug- 
gesting three-sublattice-ordered ground state in this regime. 
El 

The situation becomes more complex if t/U ~ 0(1), since 
in this regime the two-site exchange Hamiltonian is not suf- 
ficient to capture the essential physics and higher ordered 
contributions such as ring exchanges should be taken into 
account. In this paper, we focus on this regime and con- 
sider the SU(3) ring-exchange model on the triangular lattice. 
Motivated by the perturbative studies of the SU(3) Hubbard 
model at 1 /3 filling, we include the "ferromagnetic" three-site 
ring exchanges and the "anti-ferromagnetic" four-site ring ex- 
changes in the model whose Hamiltonian is given in Eq. (01. 
With the interplay between the ring exchan ges, we conjecture 
that the Quantum spin liquid (QSL) states Il6l \vh can arise 
due to the strong frustration. 

In this work, we first study the ordered phases using the 
site-factorized ansatz. fl8ll We find that the phase diagram con- 
tains the three-sublattice ordered phase, similar to the phase 
found in Ref. [l5l and the ferromagnetic state. We further 
study this model using slave-fermion trial states. After per- 
forming numerical full optimization of the trial energy, we 
consider the three ansatz states such as the uniform zero-flux 
and 7r-flux gapless (Fermi-surface) spin liquid state and the 
trimer (plaquette) state. Furthermore, we consider several 
pairing instabilities and find that there is a pairing instabil- 
ity of the zero-flux spin liquid state toward a /-wave (nodal) 
spin liquid state; there is a pairing instability of the uniform 
7r-flux spin liquid state toward a interesting exotic s-wave spin 
liquid state with two flavors of fermions paired up while one 
flavor of fermions remains gapless. 

The paper is organized as follows. In Sec.|II]we define ex- 
plicitly the model Hamiltonian we will study in the paper. In 
Sec. Ill Al we use the site-factorized ansatz to study the ordered 
states in this ring-exchange model. In Sec. Ill Bl we use the 
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slave-fermion representation to rephrase the SU(3) Hamilto- 
nian in terms of the three-flavor fermionic Hamiltonian and 
perform the fermionic mean-field treatment of the model. We 
further study the possible pairing instability in the spin liquids 
regime. In Sec. [HI] we conclude with some discussions. 



II. SU(3) MODEL WITH RING EXCHANGE TERMS 

The model Hamiltonian we consider is 

A 

+K 4 Y, [ p i234 + H.c], (3) 

with # _ # running over all the bonds on the lattice; A run- 
ning over all the triangles, up- and down-triangles, on the lat- 
tice, while (123) € P are the sites on the triangles labeled 

counterclockwise; running over all the rhombi and 

(1234) are the sites on a rhombus labeled counterclockwise. 
P\2 is the nearest-neighbor two-site exchange operator, P123 
is the three-site spin ring exchange operator, and P1234 is the 
four-site spin ring exchange operator, Fig. [T] The couplings 
J, -K3, and K4 can be obtained from perturbative analysis 
of SU(3) Hubbard model at 1/3 filling and the leading-order 
terms are 
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Without the ring exchange terms, previous studies of the site- 
factorized ansatz on the tr iang ular lattice predicted a three- 
sublattice ordered state [18-20] which was recently confirmed 
by Density Matrix Renormalization Group (DMRG) and infi- 
nite Projected Entangled-Pair States (iPEPS) analysis lfl5ll . 

Recently, Ref. |2l] did variational studies on the SU(3) 
model with three-site ring exchange, and they found the fer- 
romagnetic state and the three-sublattice ordered state in the 
regime of ferromagnetic three-site ring exchanges. On the 
anti-ferromagnetic side of the three-site ring exchanges, they 
found an interesting d x + id y spin liquid state with a gapless 
parton Fermi surface. 

However, the situation is not clear while both the three-site 
and four-site ring exchange terms are included. Here in this 
work, we will focus on the regime where the three-site ring 
exchange is ferromagnetic and the four-site ring exchange is 
anti-ferromagnetic, motivated by the perturbative studies of 
SU(3) Hubbard model. In order to study the ordered phases, 
below in Sec. Ill Al we first present our studies of ordered 
phases using the site-factorized states. Later in Sec. Ill B I we 
will present our studies using the slave-fermion trial states. 




\ e 3 = -(ei + e 2 y^ 
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FIG. 1. The triangular lattice showing the vectors §^=1,2,3 used in 
the text. The shaded region represents a single rhombi and we label 
each site counterclockwise from 1 to 4. The angle 9 between each 
vector e is 2n/3. 
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FIG. 2. The phase diagram using the site-factorized states. The red 
crosses represent the three-sublattice ordered states with each on- 
site vector, Eq. ifTJ, mutually orthogonal to each other and the closed 
blue circles represent the ferromagnetic state. The dashed line is the 
tentative boundary between the two phases. The boundary between 
the ferromagnetic state and the three-sublattice ordered state can be 
determined analytically (see text) and consistent with the numerical 
studies. We note that this model contains a four-sublattice ordered 
state for strong four-site ring exchange, roughly when K4/J > 1.4. 
In such state, the numerical site-factorized states indeed show a four- 
sublattice periodicity. However, the four-sublattice ordered state is 
hard to be characterized analytically. 



A. Site-factorized state studies 

In this subsection, we consider the site-factorized state lfl8ll 
defined as 



with 



\Xj) = aj\x)j +bj\y)j +Cj\z l3 



(7) 



(8) 
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where we fix the overall phase by setting the phase of aj to be 



zero such that a,j £ K and bj 
Ic-I 2 = 1. 



£ Z and 



We numerically calculate the site-factorized ansatz state en- 



ergy, E sf 



H 



SU(3)\S) 



on a 3 x 3, and a 6 x 6 trian- 



gular lattice and determine the optimal state in the parameter 
regime. The phase diagram is shown in Fig. [2] The red crosses 
represent the states with three-sublattice order with sublat- 
tices labeled as A, B, and C. The three-sublatticed states are 
characterized with each on-site vector (\Xj^A), \Xj^B), an d 
\Xj,zc)) mutually orthogonal to each other and the closed blue 
circles represent the ferromagnetic state. We note that this 
model shows four-sublattice ordered state with strong four- 
site ring exchange, roughly when K4/ J > 1.4. In the regime, 
the numerical values of the site-factorized states show a four- 
sublattice periodicity. However, the four-sublattice ordered 
state is hard to be characterized analytically. 1I22I l23tl 

The boundary between the three-sublattice order and the 
ferromangetic phase can be understood analytically. The 
three-sublattice state is a state with each on-site vector mu- 
tually orthogonal to each other. The state vector for the three- 
sublattice state has the form \i/j) u = \ x )j&A §5 \y)keB £3 
|z); e c with j, k, I around a triangle. Since the state vec- 
tor at each site is orthogonal to each other, the energy is 



H 



SU(3)/3-sub 



0. On the other hand, the state vectors 



of a ferromagnetic phase are aligned with each other, say 
\z) for all sites. The energy in the ferromagnetic phase is 
(Hsu {3)) fm = 3J — + 6K4. The boundary between 
these two phases can be determined analytically from the con- 
dition 3 J — 4A'3 + 6K4 = 0. For K4 = 0, the transition point 
is at 3/4, which is consistent with our numerical studies. 



The exchange operators in terms of fermions are 

p _ f »t fP fP] fa 
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where a, /3, 7, rj = x, y, z. 

The Hamiltonian, Eq. (f3]l, can be re-expressed as 

Hsu { 3) ^EE/r'/f/f/r 
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Below, we will calculate the trial energies using the slave 
fermion trial states. We start from the case without consid- 
ering pairing instability. We find that the the main competing 
states are what we call the "trimer" state, the uniform 7r-flux 
spin liquid state, and the zero-flux spin liquid state. Later, 
focusing on the regime in which spin liquid states are the op- 
timal slave-fermion states, we consider the possible pairing 
instabilities. We find that there is a possible pairing instabil- 
ity of the zero-flux spin liquid state toward a /-wave (nodal) 
spin liquid state and a pairing instability of the 7r-flux spin liq- 
uid state toward an exotic s-wave spin liquid state with two 
flavors of fermions paired up while one flavor of fermions re- 
mains gapless. 



1. Without pairing instability 



B. Slave-fermion trial states and energetics 



In this subsection, we follow the approach similar to the 
one outlined in Ref. J24I1 for the spin S — 1. We write the 
spin operators in terms of three flavors of fermionic spinons, 

r, 



s? = -*E e 

Pa 



ap-y j/3t fl 
J j J j j 
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In this section, we focus on the non-magnetic trial states. 
When we perform numerical calculations, we relax the con- 
straint of the fermion number for each flavor to be 

1 

3' 

A convenient formulation of the mean field is to consider a 
general SU(3)-rotation invariant trial Hamiltonian 



(fj fj ) trial — • 



(15) 



H 



trial 



EEh 



(jk) 



fee 



EE/o/; vi ./7 



" - /r./r • h.c. 



(16) 



with a, /3, 7 £ {x, y, z} and j is the site label. Rewriting 
the spin operator in terms of fermionic spinons enlarges the 
Hilbert space. To recover the physical subspace, a local con- 
straint on the fermions has to be enforced, 



E/f/; 



1. 



(10) 



with tjk being the hopping amplitude, 8jk being the phase 
of the hopping tjk in different mean-field ansatz states, and 
fij being the chemical potential which can be used to satisfy 
the constraint, Eq. (15[ . With the trial Hamiltonian above, 
we can find the ground state and use it as a trial wave func- 
tion for the Hamiltonian Hsu(3)< Eq- ©. After performing 
"complete" Wick contractions and ignoring the constant pure 
density terms, the trial energy can be expressed as 
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Above we defined (*? fc )* = {ff f^trial- 

The slave-fermion trial states which conserve the transla- 
tional symmetry we consider in this work are the uniform 
zero-flux state and uniform 7r-flux state. The fermions in both 
of these trial states hop isotropically on the lattice, tjk ~ t, 
and therefore the expectation values of \jk are expected to be 
isotropic, x% — X° — const. Below we list the numerical 
values of the energy per site in the two states 



jpMF 



E. 



MF 

(f> — 7T 



-0.7672 J + 0.4482A" 3 - 0.4521^4, (18) 
-0.4395 J + 0.8352/V3 - 0.7949^4- (19) 



We can observe that the optimal translationally invariant 
slave-fermion state favors = state when K% > and 
4> = n state when K3 < 0. The distinguishability between 
these two uniform-flux spin liquid states can be related to the 
origin of the /V3 which arises from the third-order perturba- 
tion of the SU(3)Hubbard model at 1/3 filling. However, K4 
does not distinguish these two uniform-flux (U(l)) spin liquid 
states from this perspective. 

Besides the translationally invariant state, we also consider 
what we call the "trimer" state. Fig. [3] shows one example 
of the configuration of such a state in which the non-zero tjk 
form non-overlapping trimer covering of the lattice. These 
states break translational invariance, and any trimer covering 
produces such a state. Such states can have lower Heisenberg 
exchange energy. The occupied bonds attain the maximal ex- 
pectation value which is found analytically \x% \max — ti° = 
1/3. Their contribution can be sufficient to produce the low- 
est total energy and such states are expected to be the lowest- 
energy states with K3 = and K4 = 0. 



E 



MF 



J + 0.5926^3 — 0.3704^4. 



(20) 



The energies of different slave-fermion trial states are func- 
tions of A'3 / J and K4/ J and it is expected that different opti- 
mal trial state is realized in different parameter regime. In or- 
der to clearly show the cross between the energies of different 



I 

mean-field ansatz states as functions of K3/ J and K4/ J, we 
plot Eqs. (Q~8]l-(|20]l with J = 1 in the limit of either K 3 = 
or K 4 = 0. 

Figure|4]shows the energies of different mean-field states as 
a function of K3/J with K4 = and as a function of K4/J 
with K3 = 0. For the former, Fig. |4(a)| clearly shows that at 
A'3, K4 = the trimer state is the lowest energy state fol- 
lowed by the uniform zero-flux state and 7r-fiux state. When 
K3 is gradually increased, the energy line of the zero-flux 
state crosses that of the trimer state and the zero-flux state be- 
comes the lowest energy state at K3/J ~ 1.61. On the other 
hand, for the latter, Fig. |4(b)| shows that the energy line of tt- 
flux state first crosses that of the zero-flux state at K4 ~ 0.96 
and then crosses that of trimer state. The zero-flux state be- 
comes the lowest energy state at K4 ~ 1.32. For a complete 
phase diagram, we numerically determine the ground states in 
the K3 — K4 parameter regimes and the result is summarized 
in the mean-field phase diagram, Fig. [5] 

In order to check if the mean-field ansatz states we consid- 
ered are sufficient to describe the physics in this model, we 
perform numerically "full optimization" of the mean-field en- 
ergy, Eq. ( TT7b , on a triangular lattice with 100 x 100 3-site unit 
cells, by treating Xjfc-s an d #jfc-s as varying variables. In the 
numerical optimization, there are totally 18 variables, 9 Xjk 
and 9 Ojk, and we take tjk = 1, fif = \x. Numerics suggest 
that the above trial states are the three optimal states. 

Before leaving this section, we want to remark that the 
trimer state is a singlet state around a triangular plaquette, and 
we can write down the exact singlet wave function in a closed 
form as 



E 

a,/3,7 



£ a/37 



|a,ft7>, 



(21) 



with a — x,y, z. With the trimer wave function, we can cal- 
culate the energy per site 



1 
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27 
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FIG. 3. Illustration of the trimer state. In slave fermion picture, the 
fermionic spinons only hop around each triangular plaquette and we 
can focus on each triangle separately. 

For if 3 = and ii 4 = 0, the exact trimer state energy is very 
close to the variational energy of U(l) spin liquid which is 
— 0.34 J|f2Hl, but much higher than the three-sublattice order 
state energy obtained by DMRG which is roughly —0.678 J. 
lfl5ll However, such a plaquette state can be stabilized with 
finite K3 before reaching the ferromagnetic state. To see 
this, we can compare the energy of the ferromagnetic state 
as shown analytically in the end of Sec. Ill Al with that of the 
exact trimer state energy. Knowing the energy of the ferro- 
magnetic state, Efm = 3 J — 4Ks + 6K4, we can see that the 
trimer state indeed has the lower energy than that of the fer- 
romagnetic state when K 3 /J < §§ + f-fCj — 0.86 + 1.5K 4 . 



(a)Energies of different trial states energies as a function of K3 while 
K 4 = 

Emf 




(b)Energies of different trial states as a function of K4 while K3 = 



2. With pairing instability 

So far, we have ignored the possible pairing instabilities in 
the Fermi surface spin liquid states discussed above. Here 
we want to address this issue. We now know that besides the 
trimer state, there are actually two spin liquid states-zero-fiux 
spin liquid state and 7r-fiux spin liquid state. Both of these 
spin liquid states are gapless and contain a single or multiple 
parton Fermi surfaces. Focusing on these regimes, we take 
the pairing mechanism into consideration and the trial Hamil- 
tonian becomes 



H, 



(Jk) a,P 



at 



ft 



a a/3 f at f /3t 
^jk J J Jk 



H.c. 



(23) 



with the constraint, n" 



(ffff) = 1/3 and t% = ±t jk 
for flux = 0/ir. For clarity, from now on we will replace 
the bond labeling (jk) by (r, r + e v ), with r running over all 
lattice sites and e v —x t 2. 3 are shown in Fig.Q] We abbreviated 
the sum over all e y= \ : 2, 3 as el 

In these regimes, the optimal state are uniform flux states 



FIG. 4. (a) Illustration of energies of different slave-fermion trial 
states as a function of A*3 with K4 = (b) Illustration of energies of 
different slave-fermion trial states as a function of K4 with A3 = 0. 
When K3 — and K4 = 0, the optimal state is the trimer state 
followed by the zero-flux spin liquid state and the 7r-flux spin liquid 
state. When K3 increases while K4 = 0, the zero-flux spin liquid 
state becomes the lowest-energy state at (K3 > 1.61, K4 = 0) as 
shown in (a). On the other hand, when K4 increases while K3 = 0, 
the energy line of 7r-flux spin liquid state first crosses the energy line 
of the zero-flux spin liquid at (A'3 = 0, A'4 ~ 0.96) and then crosses 
the energy line of trimer state to become the lowest energy state at 
(A3 = 0, K4 ~ 1.32). For general cases, the complete mean-field 
phase diagram is shown in Fig. [5] 



and the uniform expectation values of hopping functions, 
(/"t( r )/0( r + e^)) = 5 a ^x a *(^)- Furthermore, in these 
regimes, the hopping functions, x Q (e„), should be real and 
are numerically confirmed. Below, we consider two pairing 
cases: Case (1) corresponds to pairing within the same flavor 
of fermions. This requires the orbital angular momentum 
quantum number to be I = 1, 3, ... corresponding to p x + ip y , 
/-wave, ...pairing states; Case (2) corresponds to BCS-type 
pairing with different flavor of fermions. This pairing requires 



I = 0, 2 



which suggests uniform trial hopping amplitude tjk — t = 1 states. 12 ll 12411 Below, we discuss the mathematical set up for 



corresponding to s-wave, d x + id y ,... pairing 
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FIG. 5. The phase diagram of the mean-field ansatzes. The lines are 
the tentative boundaries between different phases. We note that the 
zero-flux and 7r-flux spin liquid states are both U(l) Fermi-surface 
spin liquid states. The difference is that the zero-flux spin liquid 
state only possesses a single Fermi pocket in the center while the 
7r-flux spin liquid state possesses two Fermi pockets near the hexag- 
onal Brillouin zone. The exact wave function of the trimer state can 
be written down explicitly, Eq. J21b . and we can calculate the corre- 
sponding energy exactly. 



each case separately. 

Case (1): pairing ansatz with A a/3 = 5 a @ A a . 

The pairing functions we consider are (/ a (r)/ /3 (r+e 1/ )) = 
5 a/3 A a (e u ). By symmetry arguments, the model conserves 
spatial rotational symmetry and as already explained above, 
the angular momentum of pairing functions should be odd, 
I = 1, 3, .... corresponding to p x + ip y , /-wave pair- 
ing, The ansatz can be further simplified to be x Q (e*i) = 

X Q (e 2 ) = X Q (e 3 ) = X a , and A Q (e 3 ) - A a (e 2 )e il4 = 
A a (ei)e d 28 = A a e ll2e , where 9 is 2tt/3 shown in Fig.Q] 

In addition, because the vectors / = {f x , f v , f z } and /t = 
{f x \ f , f } transform as a three-dimensional vectors un- 
der spin rotation, we expect \ x = X v = X z = X an d 
A x = A y = A z = A. In this pairing scheme, the S77(3) 
is broken down to SO(3). It is straightforward to diagonalize 
the mean-field Hamiltonian 



(24) 



a keB.Z. 



where a a (k) are Bogoliubov quasiparticles satisfying the 
transformation 



r(k) 
/-t(-k) 



a* 
k 
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'k 



a Q (k) 

3 Q t(k) 



(25) 



with 



vc = - 



1 - 



ga(k) 

E a (k) 
E a (k) 



(26) 



The ground state can be written as 

i g ^>= n n[<+</ Qt ( k )/ Qt (- k )] i^),(27) 

a.=x,y,z k 

where we define 

e Q (k) = -^2co S (k-e)-/i, (28) 

e 

A Q (k) = ^iAsin(k-e), (29) 



E a (k) = ^/(e«(k)/2) 2 + |A Q (k)p. (30) 

We can see there are three degenerate bands in this case. 

Case (2):pairing ansatz with A aa = 0, A al3 \ a ^p ^ 0. 

We consider the pairing function of the form, (f a {v)f* (r+ 
e v ))a^i3 = A Q,3 (e„). It seems there are three pairing func- 
tions we need to consider, A xy , A yz , and A zx , but in this 
SU(3) symmetric model, we can perform a global gauge trans- 
formation to make A yz , A zx = as long as the length of the 
vector formed by As is conserved, | A xy \ 2 + \A yz \ 2 + \A ZX \ 2 = 
A 2 = constant. lfll l25l - l29ll If this pairing state is energetically 
favored, there is always one gapless fermion left in this SU(3) 
system which we choose to be f z in this paper. From now on, 
we will set A yz = A zx = 0, while A xy = s/3A. In this 
gauge choice, the symmetry breaking process is more appar- 
ent. The symmetry breaking is SU(3) -> SU(2) <g> U(l). The 
SU (2) symmetry is generated by the psudo-spin doublet f x 
and f y and the U(l) is generated by the gapless f z fermion. 
After Bogoliubov transformation, the ground state is 

IGS) =nt< y + vl v r\k)f y \-k)]\vac) ® (31) 

k 

® ]T f z Hk)\vac), (32) 

where the second part, Eq. d32l , is the wave function of the 
free f z fermion, and u xy and v xy have the same expression as 
in Eq. (fJSJl with ^(k) = ^(k) = £ z (k) = ^(k) the same to 
Eq. (HU) and 



A xy (k) = ^\/3Acos(k-e), 



E xy (k) = V«(k)) 2 + |A^(k)| 2 . 



(33) 
(34) 



We can see in the SU(3)-symmetric point, the energy bands 
always show one gapless branch corresponding to one flavor 
of gapless fermions. 

With the two pairing ansatzes above, we focus on the 
regimes where the zero-flux spin liquid state and the 7r-flux 
spin liquid state are realized as shown in Fig. [5] We again 
perform full Wick contractions and ignore the pure constant 
density terms. The trial energy expression after Wick contrac- 
tions is too complex to write out explicitly. We test all differ- 
ent pairing ansatzes above and numerically calculate the trial 
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K 4 /J 



3 -- 



2 -- 



1.32 ' 



7l-flUX 

+ s-wave pairing 



zero-flux 



Trimer 



Zero-flux + f-wave pairing 



t 2 

1.61 



K 3 /J 



FIG. 6. The mean-field phase diagram with pairing instability. The 
lines are the tentative boundaries between different states. We note 
that the line separating the 7r-flux and zero-flux spin liquid states and 
the line between zero-flux and trimer state are the same as shown in 
Fig. [5] Compared with Fig. [5] the phase diagram shows that the zero- 
flux spin liquid state has a possible pairing instability toward a /- 
wave gapless (nodal) spin liquid state, and the 7r-flux spin liquid state 
has a possible pairing instability toward a s-wave spin liquid state 
with two flavors of fermions paired up while one flavor of fermions 
remain gapless. 



energies with the optimal pairing functions A in the uniform- 
flux spin liquids regime on the triangular lattice with 300 x 300 
sites. 

The result is summarized in Fig. [6] We find that roughly in 
the regime where the ferromagnetic three-site ring exchanges 
slightly dominant, the zero-flux spin liquid state (gapless spin 
liquid with parton Fermi surface) has the pairing instability to- 
ward a /-wave gapless (nodal) spin liquid. When the four-site 
ring exchange K4 is strong enough, the pairing instability can 
be suppressed and the optimal state is the zero-flux spin liq- 
uid with Fermi surfaces. However, when we keep increasing 
K4, the optimal state becomes the uniform 7r-flux spin liquid 
state. Interestingly, in the 7r-flux spin liquid state, roughly 
when the four-site ring exchanges dominant, such a 7r-flux 
spin liquid state has the pairing instability toward an exotic 
s-wave spin liquid state with f x pair with f y , which forms a 
psudo-spin singlet, while f z remains gapless. We note that 
we numerically find no pairing instability of the uniform spin 
liquid states toward p x + ip y spin liquid states or d x + id y spin 
liquid states in the focusing regime where K3, K4 > 0. 



III. DISCUSSION 

We study the SU(3) ring-exchange model with "ferro- 
magnetic" three site ring exchanges and "anti-ferromagnetic" 
four-site ring exchanges. We first use the site-factorized 
ansatz to study the model and find the three-sublattice ordered 
states, ferromagnetic states in a large regime of K3 — K4 pa- 
rameter regime. 

In the slave-fermion trial states studies, we find the main 



competing states are the trimer state, uniform 7r-flux spin liq- 
uid state and the zero-flux spin liquid. The trimer state is 
strongly suppressed by increasing the strength of the ring ex- 
changes. We also find that the zero-flux state has a possible 
pairing instability toward a /-wave gapless (nodal) spin liq- 
uid and the 7r-flux spin liquid state has a pairing instability 
toward an exotic s-wave spin liquid state with f x pairing with 
f y fermions while f z fermions remain gapless. 

These gapless spin liquid states have different properties 
and can be distinguished, at least in the mean-field picture. 
For clarity, below we will consider specifically the spin corre- 
lations, (Sj ■ Sk)conn and the (nematic) correlations of diago- 
nal elements of tranceless quandrupolar tensor, (Q°' a Q%) C onn 
with Qf = (SfSf + SfSf)/2 - 25 aP /i. Above we de- 
fined the connected correlations (O^ Ok) conn = (OjOk) — 

(Oj) (Ok)- Before jumping into the discussions of the prop- 
erties of different spin liquid states, we first note that the cor- 
relation functions of the diagonal elements of the quadrupolar 
tensor can be expressed as 



(Q" a Qk a ) 
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(fffk) 



where above we used the identity, Q c - 



(35) 



gap _ jetfc 



-i 

with the constraint 5Z a /j*Vj* = 1-JHll We can see the 
first constant is universal for different spin liquid states, but 
the second contribution is qualitatively different in different 
spin liquid states can be used for characterizing different 
spin liquid state. For simplicity of discussion, we define 
(Qj a Q^ a )-- We will also discuss the thermal properties of 
each spin liquid state. 



Properties of uniform- flux (U(l)) spin liquid states: The 
uniform zero-flux and 7r-flux spin liquid states are both U(l) 
(Fermi-surface) spin liquid states. The zero-flux spin liquid 
possesses a single parton Fermi pocket in the center and the 
7r-flux spin liquid possesses two parton Fermi pockets near 
the hexagonal Brillouin zone. In these two different U(l) spin 
liquid states, 



(Sj ■ Sk) c J m 



(Q«Qt) mS ~ 

1 + cos [(kp R - k FL ) • (r k - rj)] 



,(36) 



where the k FR y L represent the momenta of the right patch 
and the left patch of the Fermi surface for an observable 
direction. Since the two Fermi surface spin liquid states 
have different geometric information of Fermi surfaces, the 
corresponding wave vectors are quantitatively different and 
can be detected by studying the corresponding spin structure 
factors. Since the uniform-flux spin liquid states have gapless 
parton Fermi surface(s), we expect to see linear-temperature 
dependent specific heat (C v oc T) and thermal conductivity 
(ft ex T). 

Properties of f-wave gapless (nodal) spin liquid state: This 
spin liquid state possesses gapless nodal points. Interestingly, 
even though the pairing break the original SU(3) symmetry, 



8 



the SO(3) rotational symmetry related to the spin rotation is 
still preserved (because of the fact that / and p with / = 
(f x , Pi f z ) transform as three-dimensional vectors). As far 
as the low-energy physics is concerned, this spin liquid state 
possesses 



" fC, I rr\<r>t 



aaynf 



1 + cos[k c ■ (r k - rj)] 



,(37) 



with k c being the vectors connecting different Dirac points 
and we ignore all the pre-factors of each term in the above 
equation. Since this spin liquid state contains gapless nodal 
points, it should possess square-temperature specific heat 
(C v oc T 2 ). 

Properties of the exotic s-wave spin liquid state: This spin 
liquid state possesses short-ranged spin correlations due to 
the superconducting gap in the x-y channel. In order to de- 
tect the gaplessness of such a spin liquid state, we can mea- 
sure the correlation functions of different diagonal elements 
of the quadrupolar tensor. For example, the Q zz correlation 
will show the power-law behavior 



(Qf Qf) m ' 



1 +cos[(k| R -k| L ) ■ (r k -rj)] 
|r k - r;| 3 



,(38) 



with kp R y L being the momenta of the right patch and the left 
patch of the f z Fermi surface. But the correlations related to 
qxx an( j Qyy exponentially decay. As for the thermal proper- 
ties, we expect to see C v oc T and k oc T due to the gapless 
f z parton Fermi surface. 

We remark that in the slave-fermion trial states studies, we 
only focus on the non-magnetic trial states and fix the number 
density of each flavor to be equal, nf = (ff ff) = 1/3. In 
principle, we can also consider the magnetic ordered states by 
making the number density of each flavor per site different as, 
for example, (ff //) = 1, and (/?*/«) = = (fffV). Af- 
ter Gutzwiller projection of all the fermionic mean-field states 
(magnetic ordered states and non-magnetic states) l3(ill . we 
can compare the energies of different states. 

From the perspective of numerics, since there have been 
DMRG and iPEPS studies on the SU(3) Heisenberg model of 



three-flavor fermions on the triangular lattice which found the 
three-sublattice ordered state. IU5I1 We suggest possibly the 
interesting zero-flux spin liquid state be also detected in the 
DMRG and iPEPS studies on the SU(3) ring-exchange model. 

From the view of cold atom experiments. Recently, the 
cold atom experiment demonstrated a method to be able to 
add an artificial tunable gauge potential to the system. rf3~Tll 
With the tunable gauge potential, it may be possible to tune 
the sign of the three-site ring exchanges from ferromagnetic 
to anti-ferromagnetic. In that case, for strong four-site ring 
exchange, the main competing states are still the zero-flux 
and 7r-flux spin liquid states, with the trial energies similar 
to Eqs (Q~8]l-([l9]l with K 3 < 0. However, for small K±, 
the main competing slave-fermion trial states are the uniform 
7r/3-flux spin liquid state with the trial energy, E^J^^ ~ 
-0.8992 J + 0.8343^3 + 0.02443if 4 , and the tt-Aux spin liq- 
uid state. It may be interesting to explore the regime with 
anti-ferromagnetic three-site ring exchange with K4 ~ and 
it may be possible to see a phase transition between these two 
spin liquid states by manipulating the artificial gauge potential 
in experiment. 

Another interesting theoretical outlook is how the phase di- 
agram, Fig. [6] evolves if we perturb the model away from the 
SU(3)-symmetric point. In this way, the model can be con- 
nected to other models l2lll24l[32Tl which can explain the gap- 
less spin-1 spin liquids possibly realized in Ba3_NiSb20g I33I1 
or other theoretical spin-1 models. lfl9ll2^[34l - l36ll 
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